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I. INTRODUCTION 



Femtosecond lasers are known to be powerful tools for micro- and nanomachining. In 
particular, these lasers can induce periodic modulations ("Laser-Induced Periodic Surface 
Structures" or LIPSS) on the surfaces of metals, semiconductors and dielectric samples 
at relatively moderate laser fluences^"-. Furthermore, it is possible to decrease periods of 
these structures below the laser wavelength, thus rising the precision of laser nanomachin- 
ing beyond the diffraction limit 7-10 . Applications of the ripple structures are numerous. For 
instance, it is possible to colorize metals and to control over laser marking;^ 2 -. 

For further development of these applications, it is important to better understand the 
mechanisms of ripple formation by femtosecond laser pulses. The number of the observed 
structures is, however, very large making the control over the laser parameters very com- 
plicated. In general, two types of structures can be distinguished (i) resonant structures, 
where the resulting period is correlated with laser wavelength, and (ii) non-resonant struc- 
tures, which are not explicitly connected with the laser wavelength and with the coherent 
effects. Thus, the so-called Low-Spatially-Frequency LIPSS (LSFL) are belong with the 
resonant structures 8 . In addition, much larger parallel ripples of several micrometers, can 
be considered to belong with the non-resonant structures. In this case, a large number of 
pulses is required to provide a thick melted depth comparable to the structure amplitude, 
so that these ripples are rather connected with the capillary wave generation or surface 
stress 5 . The non-resonant structures were also explained by the self-organized processes, 
described by the Kuramoto-Sivashinsky equation^"—. With the increase in laser pulses 
and fluence, drop-like structures named "beads", and then conical structures called "black 
silicon" can be also obtained 41 —"-^. 

In this paper, we focus our attention at the resonant periodic structures with the period 
near the laser wavelength e.g. the LSFL. The classical theory of ripple formation proposes 
that scattering of the laser wave by surface roughness couples the laser wave with the sur- 
face modes, which interfere with the laser light, and thus lead to a periodic modulation of 
the absorbed energy 20 . This theory was recently confirmed by the numerical calculations 
based on the system of Maxwell equations for rough surface 21 . This scenario requires the 
presence of an initial roughness of a certain size. Both the laser parameters and the surface 
roughness are, however, often unknown. In addition, several laser pulses are frequently re- 
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quired to form a structure. In this case, the period of the energy deposition can be smaller 
than laser wavelength, as was explained by Bonse et al 22 by using the "Sipe-Drude" model. 
In this way, an explanation of the very narrow structures (HSFL) appearing after numerous 
pulses 23,24 was proposed. In addition, Tsibidis et al. 25 recently investigated the cumulative 
hydrodynamic effects and the corresponding surface modifications. It was found that a 
non-resonant mechanism explains the reduction of the LSFL periodicity with the increase 
of pulse number in Si. In addition, the possibility to create periodic surface modifica- 
tions with a single femtosecond pulse was demonstrated for metals and semiconductor 
materials 18 ' 22 26 . To explain the formation of the near-wavelength ripples at intense and 
reduced number of pulses, several authors have proposed the surface plasmon polaritons 
(SPP) as the mechanism responsible for the surface wave generation in semiconductors 
and dielectrics 22 - 2 -^. The surface plasmon polaritons are known to be excited on metal 
surfaces. The ripple formation mechanism for metals has been linked with the excitation 
of surface plasmon polariton by several authors 29,30 . However, the excitation conditions 
remain rather puzzling in the case of semiconductor or dielectric materials. Moreover, it 
is not clear if the SPP coupling can occur by using a single femtosecond laser pulse, since 
specific coupling conditions are required to add the missing momentum at the surface. In 
particular, gratings, snom probe, prism, defects, or roughnes s 31 ' 32 typically help to couple 
the laser wave with a surface wave mode. Thus, our study focuses on the semiconductor 
case, and the analyzed material is monocrystalline Si. 

In the case of semiconductor materials, laser-induced modification of the dielectric 
function changes material properties. As a result, the importance of a transitory metallic 
state was underlined by Bonse et al 33 . In such cases, however, the presence of a nanomet- 
ric defect (such as bubbles or nanoparticle) at the surface is required and those methods 
correspond to the well-known case of localized surface plasmon (LSP) excitation around 
isolated defects 34 . It was demonstrated, furthermore, that the transient modification of the 
solid properties follows the plasma dynamics of the free-carrier gas, due to their excitation 
by the intense laser 3 ^ 3 ^. However, a systematic study is still required for the conditions of 
SPP excitation on semiconductor surface by a femtosecond laser interaction. Under these 
excitation, electron-hole pairs are generated, but also thermal effects play a role. That 
is why a clear explanation is required for the laser parameter range, ambient environ- 
ment and sample surface conditions. To help developing the corresponding applications, 
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the resulting ripple periods should be connected with laser parameters. In this paper, we 
consider the conditions required to excite surface plasmon polaritons on semiconductor's 
surface. The developed model provides all the parameters, which lead to the SPP excita- 
tion, such as angle of incidence, laser fluence, pulse duration, and surface roughness for 
the given ambient optical properties. It is demonstrated that, under the required condi- 
tions, Si surface becomes optically active under femtosecond irradiation, and thus, SPPs 
can be excited by irradiation of a coupling device. The resulting periods are analyzed as a 
function of laser parameters. 

The paper is organized as follows. In Section [III we present the experimental protocol. 
In Section [Till we present the model and consider the modification of the optical properties 
of Si under femtosecond irradiation. In Section [iVl conditions of the excitations of the 
surface plasmon polaritons are presented. In Section |Vj, the laser parameters allowing the 
SPP excitation are examined as a function of laser fluence and laser pulse duration. Then, 
the required minimal roughness thickness is discussed. Finally, the calculated periodicities 
are compared to the experimental values as a function of laser fluence. The evolution of 
the ripple period is analyzed as a function of the angle of incidence and laser polarization. 

II. EXPERIMENTAL DETAILS 

The micromachining experiments were performed by using a Ti-sapphire laser (Hurri- 
cane model, Spectra-Physics) that was operated at 800 nm, with an energy of 500 \iJ, a 
repetition rate of 1 kHz and a laser pulse duration of 100 fs. Laser irradiation of silicon 
surface was carried out in a vacuum system with a pressure of 5 1CT 5 to 1 1CT 5 mbar. This 
low pressure considerably reduces the redeposition of unwanted debris from the laser ab- 
lation process. To get a more uniform laser energy distribution, only the center part of 
the gaussian laser beam was selected using a square mask of 2 x 2 mm 2 . A spot about 
35 x 35 nm 2 area was obtained projecting the mask image onto the sample surface with 
a lens (f = 50 mm). Laser beam was perpendicular to the sample surface. The laser 
energy delivered to the sample surface could be attenuated by coupling an analyzer and 
a polarizer and completed by a set of neutral density filters. The analyzer rotation placed 
in front of the polarizer are controlled by a computer. The engraving results are in situ 
monitored by a CCD camera. The number of pulses is controlled by triggering a Pockels 
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Physical meaning 


Notation 


Value 


Unit 


Reference 


One-photon absorption coefficient 




2wlm v / e 00 (u))/c = 1.021 10 5 


m _1 


(d, g) 


Two-photon ionization rate 




0.1 10~ 9 


m.W- 1 


(c) 


Impact ionization probability rate 


Si 


3.6 l0W e -E a /k B T e 


s- 1 


(b, e, f) 


Auger recombination rate 


C 


3.8 10" 43 


m 6 .s _1 


(b) 


Recombination delay at high density 


TO 


6 


ps 


(a) 



TABLE I: The calculation parameters for c-Si under 800 nm irradiation. References: (a) 
Ref.- 3 - 7 , (b) Ref. 38 , (c) Ref. 39 , (d) Ref. 40 , (e) RefA (f) RefA (g) RefM. 

cell, thus reducing the repetition rate of the laser pulse to 5 Hz. We irradiated a <100> 
monocrystalline silicon (c-Si) wafer by one or several laser pulses at fluences of 0.5 J /cm 2 , 
0.8 J /cm 2 , and 1.15 J /cm 2 . Two series of experiments were performed, (i) At very low 
(one or two) number of pulses, the angle of incidence has been kept normal to the surface, 
(ii) At N = 10 pulses, the angle of incidence and the laser polarization have been varied. 



III. MODELING DETAILS 



Femtosecond laser can promote carriers from the valence band of a semiconductor to 
the conduction band leading to free-carrier absorption. In our model, the number density 
of the carriers in the conduction band is calculated by solving the following equation 



dt 



- V ■ (k B T e fi e Vn e ) = G e -R e (1) 



hui 2fko 



n e +n v 



is the 



where n e denotes the electron-hole pair density, and G e 
gain of free-carriers per time unit and volume unit (m~ 3 .s _1 ). n v is the quantity of valence 
band electrons, and n e is the quantity of conduction band electrons. 

Both one-photon inteband cross-section (ai) and the two-photon cross-section (cr 2 ) are 
used in the model (Tabled]). The conduction band can be also populated due to the electron 
impact ionization (avalanche process). The corresponding coefficient 5i is also given in 
Table [|J R e = - J/ e 1 is the loss of conduction electrons by Auger recombination, where 
the recombination time t is equal to 6 ps in our calculation s 37 ' 44 . 

The initial density of free-carriers present in the conduction band is 1.84 10 9 cm -3 at a 
temperature of 300 K. k B is the Boltzmann constant. In the near-ablation regime, using 



low energy photons (1.5 eV in our case), the number of excitable electrons is limited to 
the ones available in the valence band. Even in ablation regime at the considered laser 
intensities, less than one electron per atom is usually promoted to the conduction band 44 . 
Thus, the number of the excitable valence band electrons is described by n = p Si = 
5 1 22 cm -3 , equal to the density of the Si lattice. During the excitation, the number of 
valence band electrons n v is therefore calculated by n v = n — n e . 

The free-carrier mobility is described by /i e = where v = 1.510 14 s _1 is the free- 
carrier collision frequency. Collision frequency is adjusted in agreement with melting 
fluence and melted depth given by Bonse et al 45 , consistent with Monte Carlo simulations 
of collision frequency in Si 4 ^. m e is the optical mass of electron-hole pairs, which is equal 
to 36 m e = 0.18m e0 , where m e0 is the electron mass. The density of electrons is calculated 
by using Eq. flTJ) taking into account thermal diffusion and Auger recombination. The 
hole temperature and density are considered equal to the ones of free-electrons, since 
the contribution of the electron-hole pairs to the absorption is taken into account by the 
optical mass in dielectric function. 

Laser energy absorption is calculated as follows 3 ^ 4 - 

^ = -a/,/ - {ail + a 2 I 2 )^—, (2) 

where / is the local intensity. Intensity at the surface is given by I z= o(t, x) = [1 — R (x)] I (t, x)-^-^ 
and Io(t,x) = ^v/^p-e~^te) e - ^^ 1 ) . F denotes the maximum fluence reached dur- 
ing the interaction. t = in our calculations. Spot size w and pulse duration r are 
respectively defined at the FWHM of spatial and temporal gaussian distributions. Thus, 
a x = 2 ^° n2 and a T = 2y/ ^ ln2 - The free-carrier absorption is described by aj cr = ^f 1 where 

n 2 = 0.5 ^— 3fte + \J^e (sf cr ) 2 + 3m (sfcrf^j and Sf cr = 1 — ^y+I^ describes the 
dielectric response of the free-carriers. 

During the interaction, the surface reflectivity changes depending on the angle of in- 
cidence 6 and on laser polarization. We take into account the contribution of each layer 
to the surface reflectivity by using a recurrence formula given by the transmission matrix 
models 

Ri,k 



"M+i + r i+1<k e 2i ^+ 1 



1 + r M+1 r i+lifc e 2 *+i 

where 4> i+1 = 27Th ^ +1 ^/e i+ i cos 9 i+ i is the phase shift induced by the i + 1 layer, h i+ \ is the 



thickness of the layer % + 1, and 6 i+ i is the angle of incidence at layer i + 1 given by Snell- 
Descartes law = arcsin (^J sin 0$ j . r i)i+ i is given by Fresnel equations 50 , is the 
full dielectric function given by Eq. © calculated at the layer i, and 0j is the angle of 
incidence at the layer i. The dielectric function is calculated as follows^!. 

n 

e Si (uj,n e ,u) = l + (e 00 (u)-l) ^ . p (3) 

n e + n v uj 1 (1 + z^J 

where (a;) is the dielectric constant and depends on laser wavelength 40 . n e is the free- 
carrier density, u p is the plasma pulsation defined here by to 2 = ^ff-. oj is the laser pulsa- 
tion defined by ui = 2?£. Photo-emission and thermo-emission are neglected, since Dumber 
field strongly limits the transport and thus the loss of free-carriers at the surfac e 44 ' 49 . 

Several parameters of the above model depend on temperature. Under femtosecond 
laser irrdiation, the electron-hole sub-system is excited to much higher temperatures and 
electron heat conductivity is much larger than that of lattice. To calculate temperatures of 
carriers and lattice, a two-temperature model is used. The laser-excited zone is small, so 
that ballistic effects should be also accounted for the electron sub-system. Therefore, we 
use the following ballistic-diffusive equation 51 " 53 to describe the electron-hole subsystem. 



' <>2J '" + ^ = V (D SBD VT e ) - 2" (T e - Ti ) + % 



\v dt 2 ' dt K OD1J CJ C e v B lJ ' C e 
where D S bd = jHr [v 7 ! + 4:ir 2 Kn 2 — l] is the effective free-carrier thermal diffusion 



term, based on the ratio Kn = l - M f £ - of free-carrier mean-free-path Imfp = jjy 
and size of the excited zone L ~ (2wQ : m v /ei~/c) . Specific heat capacity of free- 
carriers is taken to be equal to the classical limit C e = \k B n e . The source term Q e = 
(tkj ~E g )^ + (2twj -E g )&£- E g 6m e ~\ « + a fcr I + E g R e - \k B T e ^ describes the 
energy of free-carriers by taking into account one-photon and two-photon ionization, the 
energy loss by electron avalanche, the free-carrier heating, the release of energy due to 
Auger recombination and last term comes from the variation of the specific heat capacity 
with time, since density is strongly modified during the pulse. 

The time of free-carrier coupling to the lattice was experimentally determined by Sjodin 
et al 39 as a function of free-electron density given by r 7 = r 70 ^1 + ^\ with n th = 
6.02 10 20 cm -3 and r 70 = 240 fs. In our calculations, the coupling rate is given by 44 
^ei = — ■ Because of the slow thermal diffusion of the lattice energy, we describe the 



temperature of the lattice T S i by a classical diffusion equation, taking into account the 
energy transfered from free-carriers as follows 



Csi 



dT Sl 



V (KsiVTgi) + j ei (T e - T Si ) 



Dt 



The specific heat capacity of Si is a function of the liquid density fraction. For solid state, 
dependence with temperature is given by relations 38 C s - S i{J-m~ 3 } = 1q6 [1-978 + 3.54 10" 4 .T - 3.68T" 2 ] 
and n s -si[W/m/K] = 10 2 [1585T~ L23 ]. For liquid state, parameters are given by^~— 
Q_ Si (T) = 1.045 l0 3 pi^ S i, where pi - Si = 2520 kg/m 3 and Wi_ Si (T) = 10 2 [0.502 + 29.3 10" 5 (T - T m )]. 
Melting temperature T m depends on the free-electron density as described by Eq. © . Dur- 
ing the phase transition, both Si heat capacity and conductivity are calculated using the 
fraction of liquid -q, and are respectively defined by Csi(T) — (1 — rj) C s si{T) + r]Ci-si(T) 
and K S i(T) = (1 — 77) k s ^ S i(T) + rjKi- S i(T). Melting is considered by using melting en- 
thalpy AH m at the melting temperature T m . The resulting thermal energy is given by 
U = J^ m C s _si(T')dT' + AH m + ff m C^si{T')dT' where U is the internal energy of the 
lattice, T is the Si temperature, T = 300 K is the initial temperature of the system, and 
AH m = 4 10 9 J/m 3 is the melting enthalpy of Si. 

Previously, two different types of phase transitions were shown to take place for Si—"—. 
The first one is thermal melting and is described by a thermal criterion based on the re- 
quired energy E m = k B T m p S i+AH m . The second phase transition mechanism is a so-called 
"non-thermal melting" due to the lattice decomposition due to a large number of carriers 
in the conduction band. The contribution of the non-thermal melting is also taken into 
account by the decrease of the band gap energy as a function of free-carrier density 38 (lim- 
ited to positive or null values) expressed by E g (T,n e ) = 1.17-4.73 10~ 4 y^3g -1.5 10" 10 ne /3 
and by a lowering of the melting temperature described by the relation 60 



m 



m 



.0 




(4) 



where = 1687 K. 



Boundary conditions for transport equations are set so that free-carriers do not leave 
the sample. The sample is 250 /xm thick, and the optical transmission has been checked to 
be zero trough the sample. 
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IV. SURFACE PLASMON POLARITONS 



Light can be coupled from free space into the surface plasmon polaritons (SPP) only by 
matching the momentum of the SPPs. This can be done via index matching 31 , or grating 
couplin g 61 ' 62 . In addition, other cases can be considered. A non-resonant excitation can 
be performed by scattering of the laser wave on surface defects or a surface roughness. In 
such a case, laser wave is scattered on a broad angular distribution, and a part of the laser 
energy couples with surface modes. Laser wave can also interact with near- wavelength 
structures as described by Mie scattering, which leads to the excitation of localized surface 
plasmons (LSP)^ 3 -. 

In each case, the excitation of surface waves requires several resonance conditions 31 * 3 ^ 4 -. 
In this part, we present theoretical conditions allowing the excitation of the Surface Plas- 
mon Polaritons (SPPs) at the laser-irradiated surface of Si. 

The excitation conditions of surface plasmon polaritons at a flat surface is that the cor- 
responding curves cross in the dispersion diagra m 31 ' 64 . The dispersion relation is obtained 
from the boundary conditions of the electric and magnetic field at the interface. The conti- 
nuity of the electric field at the interface results to the expression |^ = — g-, where £i j2 are 
the dielectric constants on both sides of the interface, and fc lj2 are the respecting momenta 
of the both sides of the interface. In the general case, this expression can be verified only 
if 

^e{e 1 )^e(e 2 ) <0 (5) 

which corresponds to a metal-dielectric interface. The continuity of the magnetic field at 
the interface leads to the expression of the SPP wave-number 



c]J e 1 +e 2 

and to a second SPP excitation condition given by 

9£e (e a ) < -sfte ( £l ) (7) 

By introducing Eq. d3]) into Eq. ([6]) and considering an interface with vacuum (ei = 1), 
we calculate the dispersion relation as a function of the free-carrier density. Figure Q] 
presents the dispersion relation of Surface Plasmon Polariton and laser wave in the case 
of Si, as a function of the laser frequency, normalized to plasma frequency. In this cal- 
culation, the band structure of Si has been taken into account by replacing £ ro (u) in Eq. 



10 19 cm 3 L jg ht |i ne 

W 2l cm~ 3 Plasma frequency 

5 10 21 cm,-'' 

50 
40 
30 

-— > 20 
W 


-10 
-20 

0.1 1 10 

u/uo p 

FIG. 1: Real part of the dielectric function of Si s(cu) = (^) 2 asa function of the photon 
energy compared to the plasma frequency uj p = J ^f^. If the plasma frequency is located 
on the left of the resonance peak, surface mode is not bounded to the surface and SPP 
are not allowed. If the plasma frequency is located on the right of the surface plasmon 
polariton resonance peak, the SPPs are allowed. 

d3]) following the measurements of Palik 40 . The wave-vector k is squared and normal- 
ized to the laser pulsation to obtain the dielectric function, and is presented as a function 
of the pulsation uj. The dispersion curve of the SPP reveals a plasmon resonance of the 
free-carrier plasma if the plasma frequency is comparable with the laser frequency e.g. if 
the free-carrier density is high. In the case of a doped semiconductor near the eigen fre- 
quency, the intersection of the dispersion curves is possible 31 . In the case of femtosecond 
laser interaction, similarly, laser-induced ionization provides free-carrier contribution to 
the dielectric function, which is equivalent to a transient doping. Therefore, the intersec- 
tion also becomes possible. The bounding of the coupled wave to the surface is described 
by the condition uj < ^| in a perfect conductor 6 -. In the case of damped material, how- 
ever, the leaky part of the dispersion relation between uj spp = ^| and uj p is allowed 64 . The 
condition on pulsation in a damped material is thus described by 

uj < uj p . (8) 

The minimal density leading to satisfy this later condition is the critical density at which 

9f?e (e) = 0. 

The coupling of far field laser wave with the SPP is also possible in the case of a surface 
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with defects or roughness 31 . In this case, the pseudo-grating period has a large thickness 
8k around its average value, and leads to a small but non-zero coupling efficiency. Be- 
cause of the laser irradiation, a coupling of few percents is sufficient to obtain a periodic 
modulation of the deposited energy by interference between laser and SPP waves 28 . 

Finally, the conditions to satisfy (Eq. [5l Eq. and Eq. [8]) allowing the excitation of the 
SPP at the vacuum- Si interface can be combined into the condition 

3fJe (e) < -1 (9) 

This criterion corresponds to a minimal free-carrier density given by N e = 4.27 10 21 
cm 3 in the case of a laser irradiation of Si using wavelength of A = 800 nm. Then, we 
determine laser parameters for which this condition is justified. 

V. RESULTS AND DISCUSSION 

A. Conditions for SPP excitation on Si 

We now demonstrate that excitation of Surface Plasmon Polariton is possible on Si un- 
der femtosecond laser irradiation, and calculate the laser parameters leading to satisfy the 
surface plasmon polariton conditions. Figure [2] (a) shows the maximum of the free-carrier 
density as a function of the laser fluence. At short pulse duration, the maximum free- 
carrier density is non-linear as a function of fluence. This effect is due to the modification 
of absorption near the critical density n cr = 3.98 10 21 cm -3 . At low laser intensity, the 
maximum density increases linearly, since the absorption is linear in this regime. Below 
critical density, one observe that the efficiency of the absorption increases, which is due 
to the significant contribution of the multi-photonic absorption. Above the critical density, 
absorption becomes limited by the surface reflectivity. Above this limit, the number of free 
carriers still increases due to high temperature of the free-carriers, leading to an interplay 
between diffusive transport and impact ionization. 

Figure [2] (b) shows that the resonance condition is not met during the irradiation with a 
long pulse duration at the considered fluences, since the quantity of free-carriers is limited 
by the low intensity. It is shown that the real dielectric function decreases linearly if laser 
fluence is near the modification threshold ( 0.2 J /cm 2 , see Ref s 45 ' 66 ' 67 ). We observe that 
the critical density is not reached in this fluence regime. The difference with Ref 36 is 
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FIG. 2: (a) Minimum of the real part of the dielectric function as a function of laser 
fluence. Blue line indicates the critical density, (b) Maximum of the density reached 
during the interaction as a function of laser fluence. Blue line indicates while 
3fte (e) = —1. Several laser pulse durations are represented. Laser wavelength is 800 nm. 



explained by the different collision time, the two-photon cross section which is 10 times 
lower here (see Ref 6 ®), and the impact ionization that we took into account. In the case 
of 100 fs pulse duration, the condition given by Eq. ([9]) is satisfied above laser fluence of 
0.7 J/cm 2 . A threshold for the excitation of SPP is then identified for a laser fluence of 
0.7 J/cm 2 , a pulse duration of r = 100 fs, and a laser wavelength of A = 800 nm. It is also 
shown that under shorter laser interaction, the condition for surface plasmon polariton 
excited is satisfied from lower fluences e.g. at 0.5 J/cm 2 if r = 50 fs and A = 800 nm. 
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From those results, a threshold fluence for SPP resonance can be defined for each pulse 
duration, above which the SPP resonance conditions are met. This result explains why 
a high fluence is necessary to induce the formation of periodic structures in single pulse 
experiment s 4 ' 26 ' 66 , since it results in a sufficient quantity of free-carriers to excite surface 
waves at the surface of Silicon. 

In Figure [3l we demonstrate the fluence threshold for the SPP resonance as a function 
of the laser pulse duration. The corresponding intensity, at which resonance occurs for 
100 fs pulse, is shown by the dashed curve. By comparison of the curves, we observe that 
the required intensity for the SPP resonance increases with the decay in the laser pulse 
duration. This effect is due to the screening and large density gradient at the surface 
resulting into a strong diffusive transport. 

Next we calculate the lifetime and the depth of the optically active zone, e.g. the 
distance under the surface where the sufficient number of free carriers are excited. Figure 
[4] shows the real part of the dielectric function as a function of time and depth. It is shown 
that under 50 fs pulse duration, at a fluence of 0.62 J /cm 2 , the excited zone is nearly 20 
nm deep and the SPPs are excited during a picosecond, which is greater than the pulse 
duration, thus leading to the excitation of SPPs in a shorter timescale than necessary for 
surface melting. Moreover, the damping length of the SPP is given by the relatior^^ 4 - 
L S pp = [2Qm (/3)]~ . The value of the damping length is contained between 500 nm and 
2 /xm if the SPP resonance conditions are met. Then, the excited SPPs propagate through 
several micrometers and can lead to periodic modulation if laser interferes with the SPPs, 
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FIG. 4: Distribution of the dielectric function as a function of depth and time, r = 50 fs, 

F = 0.62 J/cm 2 . 

as experimentally observed around defects^^2^2iZ£. 

As underlined in the previous section, the phase-matching is possible at the surface 
of Si by using a scattering configuration with a defect or a roughness. We separate the 
following cases (i) the case of the roughness, for which the size of the scattering center is 
very small compared to the laser wavelength, and (ii) the case of a defect, for which the 
size is comparable to the laser wavelength. Both situations can lead to the excitation of 
Surface Plasmon Polariton if the conditions on the dielectric function are justified. 

Figure [5] shows the distribution of the time-averaged transmitted field amplitude and 
normalized by the incident laser field intensity at the bottom of the selvedge region. This 
result has been calculated using the FDTD simulation package Lumerical 71 by irradiating 
a randomly generated rough Si surface using several control parameters: FWHM of the 
amplitude, and distance between the scattering centers. In this calculation, the roughness 
amplitude is distributed as a gaussian function between and 15 nm. The distance be- 
tween scattering centers is taken equal to 100 nm so that the scattered waves interfere 
together. The Si dielectric constant is taken equal to the value under 800 nm wavelength 
laser irradiation, in the case of low laser excitation. One observes that the amplitude of 
the field transmitted below the roughness is modulated if roughness amplitude is greater 
than 8 nm. Such a roughness is formed after a single laser pulse 7 ^ at 0.5 J/cm 2 . Thus, the 
amplitude allowing the coupling of surface waves with laser is 8 nm, which explains why 
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FIG. 5: Transmitted field amplitude below the selvedge region of Si. Roughness 
amplitude is (a) 5 = nm, (b) 5 = 8 nm, (c) 5 = 15 nm at wavelength 800 nm. In these 

simulations, 6 = 6^ = 13.64 + 0.048i. 

stricly parallel ripples are observed after two pulses or more. Conversely, the formation 
of single pulse periodic structures is due to scattering on a near-wavelength defect, which 
leads to the excitation of Localized Surface Plasmon Polaritons distributed around scatter- 
ing centers as observed by several author s 22 ' 26 ' 6 ?. High fluence single pulse experiments 
leads to the formation of concentric structures rather oriented in the direction perpendic- 
ular to the laser polarization. 
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FIG. 6: Wavelength normalized period of SPP as a function of free-carrier density at 
vacuum - Si interface when the conditions of resonance are met. 

In this section, we have theoretically demonstrated that the excitation of Surface Plas- 
mon Polaritons occurs on Si irradiated by femtosecond lasers. The excitation conditions 
are satisfied during the laser pulse if the laser intensity is high. At 50 fs pulse duration, a 
layer of about 20 nm becomes optically active and has a lifetime longer than the the pulse 
duration. We turn now to the study of the period of the SPPs excited during ultrashort 
laser pulse on Si. 

B. Effect of the experimental parameters on SPP period 

The SPP period dependency on laser intensity depends on the free-carrier density as 
follows 



where e\ and e-iiuS) are respectively the dielectric functions of the media at both sides of 
the vacuum - Si interface. By substituting e 2 (oo) with Eq. O), the periodicity of the SPP as 
a function of free-carrier density is calculated. 

Figure [6] demonstrates the period of the SPPs at the vacuum - Si interface as a function 
of the free-carrier density. The values are presented for the free carrier number densities 
required for the SPP excitation. The resulting period varies considerably with the carrier 
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FIG. 7: Experimental measurements of the LSFL periods, as a function of laser fluence. 
r = 100 fs, A = 800 nm, 6 = 0°. The periods resulting from theoretical investigations are 
also represented. Bonse (2010) refers to Ref. 4 and Bonse (2002) to ReiM 



density. The period of the SPPs is contained between 0.7 A and A, which correlates with 
the generally observed LSFL periodicitie s 4 ' 22 ' 66 . A quantitative study of the variation of 
the SPP periodicity with laser fluence is now presented, and compared to the LSFL ripples 
formed using a very low number of laser pulses. 

Figure [7] shows both theoretical and experimental periodicities. The period of the LSFL 
structures is presented as a function of laser fluence, for 100 fs pulse duration. When SPP 
resonance conditions are satisfied, the resulting SPP period tends to the laser wavelength 
when increasing laser fluence. In the optically active range (fluence is greater than 0.7 
J /cm 2 and pulse duration r = 100 fs), the calculation results agree with the presented 
experimental measurements taken from Ref s 4 ' 66 ' 72 at very low number of pulses. The 
single pulse case is explained by excitation of SPP via coupling with a surface defect. The 
case iV = 2 is explained by coupling with roughness. This result shows that the periodicity 
of Surface Plasmon Polaritons well describes the evolution of the structure period as a 
function of the laser fluence at reduced pulse number. 

In the case of 10 laser pulses, LSFL ripples are well-developed, and scattering by a 
grating model well explains the observed periodicities, which confirms Ref 73 . Figure [8] 
shows the comparison between theoretical variation of period with angle of incidence with 
experiments made using 10 pulses at various fluence s 70 ' 72 . Both directions of polarization 
are presented. The theory of scattering by a periodically structured surface leads to a 
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FIG. 8: Ripple periodicity as a function of angle of incidence and laser polarization after 
10 laser pulses. Pulse duration is 100 fs, and laser wavelength is 800 nm. 



periodicity given by 62 



where n . , 
For S polarization, 
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for polarization parallel to the plane of incidence. 
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rj + sin 9 

The variations of both measured and calculated ripple periods are explained by the varia- 
tion of rj near the critical density. Actually rj ~ 1 if Re (e) <C — 1, and rj varies between 0.5 
and 1.5 if the condition (Eq. [9]) is satisfied. This section demonstrated that the theory of 
SPP excitation on gratings agrees with the experiments and explains the variation of the 
ripple period with the angle of incidence and with laser polarization. We underline that 
the threshold fluence allowing SPP resonance on the surface of a grating is decreased with 
respect to the results presented in Figure [2] (a), since the energy absorption is enhanced 
in the presence of the grating 61 . 



VI. CONCLUSIONS 



The possibility of the surface plasmon polariton excitation on Si surface irradiated by a 
femtosecond laser pulse has been theoretically demonstrated. A sufficient number of free- 
carriers is excited from the valence band to the conduction band during the laser pulse, 
thus satisfying the SPP excitation conditions. The required ranges of laser fluences and 
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pulse durations have been identified to satisfy the SPP excitation conditions. In particular, 
SPPs can be excited by using a femtosecond laser with 800 nm wavelength, 100 fs pulse 
duration and with laser fluences larger than 0.7 J /cm 2 . The shorter temporal laser pulse 
width leads to the increase in the intensity threshold for the SPP excitation. This effect is 
due to the thermal diffusion of free-carriers and to the surface screening. As a result, a 
layer is excited, with a lifetime longer than the pulse duration and with a depth of several 
tens of nanometers. 

Furthermore, a comparison of the calculated SPP periodicities and experimentally mea- 
sured ripple periodicities allows us to conclude that the formation of periodic structures 
with a reduced number of laser pulses is due to the excitation of SPPs at the Si surface. 

The presence of a surface roughness with ^ < A leads to the coupling of the laser 
wave with the roughness. We have found that the required roughness amplitude allowing 
the coupling of laser wave with the surface is 8 nm. It is also possible to obtain periodic 
structures by scattering on defects (5 ~ A) that are present at the surface by excitation 
of localized surface plasmon polaritons. These results underline the importance of the 
surface quality in the SPP excitation and thus to the LSFL ripple formation. 

As a result of the performed analysis, the possibilities of control over the period of the 
LSFL ripples can be deduced in the regime of low number of pulses. The period can be 
reduced down to 40% by increasing the angle of incidence for S polarization, and can be 
increased up to 37% by increasing the angle of incidence for P polarization. The period 
of the LSFL ripples can be increased by about 10% by increasing the laser fluence up to 
5 J /cm 2 . Finally, one should mention that LSFL period reduction due with pulse number 
also allows to decrease the periodicity of 50%, but the corresponding mechanism is still 
under discussions^^. 
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